Spin waves in the (tt, 0) magnetically ordered iron chalcogenide Fei.osTe 
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We use neutron scattering to show that spin waves in the iron chalcogenide Fei.osTe display novel 
dispersion clearly different from both the first principle density functional calculations and recent 
observations in the related iron pnictide CaFe2As2. By fitting to a Heisenberg Hamiltonian, we 
find that although the nearest-neighbor exchange couplings in the two systems are quite different, 
their next-nearest-neighbor (nnn) couplings are similar. This suggests that superconductivity in the 
pnictides and chalcogenides share a common magnetic origin that is intimately associated with the 
nnn magnetic coupling between the irons. 

PACS numbers: 74.70.Xa, 78.70.Nx, 75.30.Ds 



All parent compounds of cuprate superconductors are 
antiferromagnetic (AF) Mott insulators characterized by 
the same local moment Heisenberg Hamiltonian [T]. For 
this reason, it is believed that magnetism is important 
for the high-T c superconductivity [2]. The iron-based 
superconductors [3j 0] share many features in common 
with the cuprates, which leads many to conjecture that 
the magnetism present in these compounds is vital for 
the presence of superconductivity. The iron-based super- 
conductors can be divided into two chemical classes, the 
iron pnictides such as CaFe2As2 and iron chalcogenides 
Fei+j,Te. Many properties of the pnictides and chalco- 
genides are similar, including similar band-structure [5] 
and magnetic excitations in the superconducting com- 
positions 6 12 . Furthermore, the magnetism in the 
pnictide parent CaFe2As2 [Fig. [TJb)] is consistent with 
first principle density functional calculations [13] , How- 
ever, the parent compound [141 115] of the iron chalco- 
genides, Fei +y Te, possesses a different AF order [Fig. 
[it a)]. Therefore, it is important to determine if mag- 
netism in these two systems can be described by a sim- 
ilar Hamiltonian. If the magnetic description between 
systems is entirely dissimilar, then it presents a serious 
challenge to many theories [T6UI9] where superconduc- 
tivity has a magnetic origin. 

By studying the spin-waves in Fei.osTe, we compare 
the magnetic couplings within the pnictide and chalco- 
genide systems. We show that although the nearest 
neighbor (nn) couplings in the two systems are very dif- 
ferent, the effective next nearest couplings (nnn) Ji are 
very similar. While our results are consistent with the 
theoretical idea that J2 is important for superconductiv- 
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ity |18j . the isotropic J2 we find in Fei.osTe is very dif- 
ferent from the anisotropic J2 yielded from density func- 
tional calculations [20 . Our results suggest that while 
the nn coupling may change, it is the nnn coupling that 
persists between different iron superconductors. 

We have used time-of-flight inelastic neutron spec- 
troscopy to determine the dispersion of spin-wave ex- 
citations in Fei.osTe (with AF ordering temperature 
T N = 68 K, see Fig. [TJd) and ref. &T$, the x = 
(non-superconducting) member of the isovalently sub- 
stituted Fei+yTei^Sez iron chalcogenide superconduc- 
tors [53]. By measuring spin- wave excitations in 
Fei.osTe throughout the Brillouin zone (BZ), we have 
used a Heisenberg Hamiltonian to determine the effec- 
tive exchange couplings of the system. Our neutron scat- 
tering experiments were carried out on the HB-1 triple- 
axis spectrometer at High-Flux-Isotope-Reactor and on 
the ARCS chopper spectrometer at Spallation-Neutron- 
Source, Oak Ridge National Laboratory, USA. We also 
used MAPS chopper spectrometer at ISIS, Rutherford- 
Appleton Laboratory, UK. For the experiment, we have 
co-aligned 6 grams of single crystals of Fei.osTe. All data 
were collected at around 10 K (<C T N ) with incident neu- 
tron energies E t = 55, 90, 180, 350, 500 and 580 meV 
with the c-axis aligned along the incident beam direc- 
tion. Since the spin-wave excitations have weak c-axis 
coupling, we integrate the excitations along the c-axis 
direction, and focus on spin waves in the (h, k) plane. 

For Fei+j,Te with modest excess iron content y, the 
magnetic structure is shown in Fig. [IJa) [HI \TE\ , which 
can be viewed as two AF sub-lattices as shown by darker 
and lighter colored atoms. We define the nn (J\ a , Jib), 
the nnn (J2a, J2b), and the next-next-nearest neighbor 
(J3) exchange interactions as shown in Fig. [It a) |20) . 
The nn magnetic exchange couplings {Ji a , Jib) are de- 
fined similarly to those of iron pnictides [Fig. [ljb)] . How- 
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FIG. 1: (a) Schematic of in-plane Fe spins displaying magnetic 
order in Fei+^Te with small y [141115] , and showing definition 
used for exchange energies, (b) Schematic of in-plane mag- 
netic order in CaFe2As2 |24| with exchange energy definitions, 
(c) Schematic showing wave vector dependence of intensity at 
various energies (for raw data see Fig. [5|. Dashed line shows 
one BZ. (d) Temperature dependence of elastic scattering at 
magnetic Bragg peak for the Fei.osTe sample. 



FIG. 2: Constant energy slices of the spin-waves as a function 
of increasing energy at 10 K for Fei.osTe. All data are normal- 
ized to absolute units with a vanadium standard, (a)-(c) col- 
lected with incident neutron energy Ei = 90 meV on ARCS, 
(d)-(f) E, = 350 meV on MAPS, (g)-(h) E, = 500 meV on 
MAPS. The dashed line in (a) shows a crystallographic BZ. 



ever, the nnn couplings (J2 a , J2b) in chalcogenid.es are 
directionally dependent as shown in Fig. [lja). 

Our Fei.osTe samples were grown using Bridgman 
technique as described before [5T]. Fei+yTei-^Sej; is 
tetragonal at high temperature and becomes orthorhom- 
bic or monoclinic (depending on x, HH [HI |22l E3] ) be- 
low Tjy. The a&-plane lattice parameters for the vari- 
ous phases remain very similar, and on cooling into the 
low symmetry phase the sample becomes twinned. We 
therefore measure the wave vector in tetragonal (h,k,l) 
reciprocal lattice units, with in-plane lattice parameters 
a = b =3.80 A, and the out-of-plane c = 6.23 A. In 
this notation, magnetic order in powder Fei+ y Te has 
been found at (0.5,0,0.5) for small y, and increasing 
y will lead to incommensurate magnetic order [UJ [TS] . 
In the present single crystalline samples, the magnetic 
order was found to be centered very close to the com- 
mensurate position at (0.485,0,0.5) r.l.u and y =0.05 
was measured with inductively coupled plasma analysis 
|21j . However, we also observed a weaker magnetic peak 
at (0.37, 0, 0.5) r.l.u attributed to a small portion of the 
sample with slightly different y. Figure [ljd) shows the 
temperature dependence of the magnetic Bragg intensity 
at Q = (0.485, 0, 0.5) r.l.u confirming T N = 68 K. 

The magnetic excitations probed by neutron scattering 
in our Fei.osTe sample are summarized by representative 
constant energy slices in Fig. [2[ The data have been 
normalized to a vanadium standard and plotted in ab- 
solute units, without correction for the magnetic form 
factor, causing the signal intensity to decrease with in- 
creased Q. Each Ei probes a different out-of-plane wave 
vector for each energy transfer, and it was found that 
data from different E^s were consistent, implying little 



i-dependence of the data over the energy range probed. 

Spin waves in most materials tend to display a mag- 
netic response centered on the magnetic Bragg position 
up to the highest energies, with successively larger rings 
with increased energy. However, we discuss below how 
the center of the excitations switch from the (0.5, 0) low 
energy position to integer positions at higher energy, 
which we interpret as the outcome of the interaction of 
competing ferromagnetic and AF exchange energies. 

At our lowest energy, 7.5 meV [Fig. [^a)], magnetic 
excitations emerge from the AF Bragg position (0.5,0) 
and other half-integer reciprocal lattice vectors [in an 
untwinned sample, magnetic peaks would not appear at 
(0, 0.5), but twinning leads to an equal intensity domain 
rotated by 90° in-plane]. As the energy is increased, the 
response spreads out in Q as expected for spin-waves 
[Figs. [2jb) andpjft]. As the energy is raised to around 
60 meV [Fig. (2Td)], there are no longer peaks at half- 
integer positions, but instead there are rings of radii 
~0.5 r.l.u which are centered on integer reciprocal lat- 
tice points. These rings are even clearer when the data 
are corrected for the magnetic form factor drop-off at 
high wave vector (see supplementary material). As en- 
ergy is increased, the radii of rings around (1,1) expand 
and those around (1,0) contract [Figs. pje)-l2jf)]. Even 
at 115 meV a ring can be seen around (1 , 0)7 which by 
225 meV contracts into a peak at (1,0) [Fig. |2jg)] before 
the disappearance of all intensity at higher energies [Fig. 
2(h)]. Corresponding cuts along the (h,0) trajectory are 
shown in Fig. [3} A schematic of the dispersion of the 
magnetic response is shown in Fig. [ljc). The data above 
100 meV in Fei.osTe have similarities to the highest en- 
ergy spin excitations observed in FeTei-^Se^ with x = 
0.27, 0.49 [H]. 
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FIG. 3: Constant energy cuts along the (h, 0) trajectory, each 
from a slice in Fig. [2] Solid lines are fits to Gaussians. 



In order to extract effective exchange energies, we fit 
spin-wave data using a Heisenberg Hamiltonian (see sup- 
plementary material for the model Hamiltonian) with 
commensurate (0.5, 0, 0.5) AF [55]. In order to yield this 
commensurate AF, there are constraints on the bounds 
of each of the magnetic exchange energies [26] . Because 
of the twinned nature of the sample, the model used is 
the sum of two equal sized domains rotated by 90° . 

To determine the dispersion curves for spin waves, the 
slices in Fig. [2]were cut along the (h, 0) and (1, k) direc- 
tions. By fitting Gaussians to many (h, 0) cuts of differ- 
ent energies like those in Fig. [3] we obtain the dispersion 
plot in Fig. |4ja) using the fitted peak positions. Simi- 
larly, (1, k) cuts were fitted to create Fig. |4|b). These two 
dispersion plots were simultaneously fitted to the disper- 
sion of the model [25], yielding the fit displayed in Figs. 
[ija)-(b). Similar conclusions about the dispersion could 
be reached by viewing the data in terms of constant- 
Q cuts instead of cuts at constant energy, but this was 
not found to be as effective for quantitative analysis (see 
supplementary material) . In Fig. |lj the intensity of the 
excitations of the model is proportional to the radius of 
the marker (which is saturated at the lowest energies to 
maintain figure clarity), to highlight the bands with neg- 
ligible intensity (also see the supplementary material for 
a zoom into the low energy part of the plots) . The pres- 
ence of almost non-dispersive bands around 250 meV is 
not clear in the Q-cuts, possibly because of averaging-out 
in Q as the bandwidth is comparable to the instrument 
resolution (along with poorer statistics at high energies). 
It is also not clear if these bands can be seen in constant - 
Q analysis (see supplementary material). 

In the fit lines displayed in Fig. |4|a)-(b), J 2 b was fixed 
equal to J2a, after it was found that these two parameters 
had very similar values when allowed to vary (see supple- 
mentary material for fit with Ji\, not fixed to j r 2a )- This 
four parameter fit leads to exchange energies of J la = 
-17.5 ± 5.7, J u = -51.0 ± 3.4, J 2 = J 2a = J 2b = 21.7 
± 3.5, J3 = 6.8 ± 2.8 meV (assuming 5=1) and fits 
the dispersion in these directions well. By further fix- 
ing J 3 = 0, the model can successfully fit the data up 
to ~ 100 meV, but the maximum band energy, CL> max , is 
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FIG. 4: (a)-(b) Solid black markers are dispersion data found 
from fitting Gaussians to form factor corrected data at many 
energies for the (h,0) and (1, k) directions respectively. Gray 
open circles (with radius indicating intensity) show best fit 
dispersion curves with fitting parameters given in the main 
text, (c)-(d) Data as in (a)-(b), but with dispersion curves 
simulated using exchange constants predicted by density func- 
tional calculations, which clearly do not agree with the data. 

underestimated by around 50 meV (see supplementary 
material for fits where J3 is fixed to zero). 

Using the fit parameters listed above, we show in Fig. 
[5] constant energy slices calculated from the resolution- 
convolved model. Here we have also considered the out- 
of-plane (c-axis) exchange coupling J z and found that 
J z = 1 meV best fits the spin-wave intensities, although 
the simulation slices otherwise do not change significantly 
with J z . The overall features of the model fit arc: (i) 
below ~30 meV, intensity is located around (0.5,0); (ii) 
at intermediate energy there are rings around (1,1) that 
grow with increasing energy; (iii) above ~150 meV the 
intensity ends in a peak at (1, 0). The data are consistent 
with the model, though the intermediate energy features 
are more grid-like than the more rounded data. 

Our fits and simulations show highly anisotropic in- 
plane nn exchange couplings with \Jib\ ^> \Ji a \i and a 
nnn exchange that is AF (energy ~20 meV) and isotropic 
J 2 = J 2a ~ Jib- The w m ax observed is between 200- 
250 meV. Comparing our results to similar high energy 
measurements of CaFe2As2 [27], which has J\ a — 50± 10, 
Jib = -5.7± 5, J 2 = 19± 3 meV and w max «200 meV, it 
is clear that the w max and values of J 2 are similar, as well 
as the presence of anisotropy in J\ in both cases plus no 
anisotropy in J 2 in either case. However, the dominating 
Ji exchange constants are -50 meV (Jib) and +50 meV 
(Ji a ) for Fei.osTe and CaFe 2 As 2 , respectively. 

Our results shed new light on the nature of the mag- 
netic state in the iron chalcogenides and its relationship 
to superconductivity. The isotropic J2 suggests that 
this nnn exchange coupling originates from the super- 
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FIG. 5: Resolution convolved simulation (using Tobyfit [28] ) 
of the Heisenberg model using the best fit parameters in the 
text plus an out of plane coupling of J z = 1 meV. Each slice 
corresponds to a slice in Fig. [2] The model has been given 
a line-width of 10 meV before resolution convolution, though 
adjusting the line-width does not make a substantial differ- 
ence. All slices are on the same intensity color scale as Fig. 
[2j with an overall intensity scale that was chosen so that in- 
termediate simulation slices had a similar intensity to the in- 
termediate raw data slices. 



exchange mechanism, and is insensitive to the lattice dis- 
tortion and variation in the ci-orbital components. The- 
oretically, it has been shown that the nnn |18j magnetic 
coupling can cause an s ± -wave pairing that induces a 
neutron spin resonance at wave vector (0.5,0.5) [2T)1 130j . 
Similar isotropic AF J2 values in iron-pnictides and iron 
chalcogenides therefore naturally explain the experimen- 
tally observed neutron spin resonance within both classes 
of iron-based superconductors 6 -12]. First principles 
density functional calculations [20] on Fei.o6sTe predict 
highly anisotropic nnn exchange interactions which are 
not consistent with our data [see Figs. 4(c) and 4(d) 
for dispersion and the supplementary material for simu- 
lation slices] , perhaps due to the complex nature of the 
orbital ordering [3TJ [32] or itinerant magnetism [33] in 
this material. 

In summary, we have shown that spin-wave excita- 
tions in the iron chalcogenide Fei.osTe can be modeled 
by a Heisenberg Hamiltonian with anisotropic (domi- 
nantly) ferromagnetic nn and isotropic AF nnn exchange 
couplings. While the nn couplings for Fei.osTe and 
CaFe 2 As2 [27] are different, we find that the AF nnn 
exchange couplings in these two classes of materials are 
not only similar in magnitude but also directionally inde- 
pendent, even though they have different AF and crystal 
structures [TU [T5] [23] . Our findings suggest that super- 
conductivity in both classes of iron-based superconduc- 
tors shares a common magnetic origin that is intimately 
associated with the AF nnn exchange couplings [18]. 
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FIG. SI: Constant energy slices of data (from Fig. 2) con- 
verted into units of magnetic response x (<1> W ) (after sub- 
tracting constant backgrounds). After the form factor cor- 
rection, the magnetic signal no longer decreases in intensity 
at high wave vector transfer, and so the magnetic pattern 
is somewhat clearer than in the raw data. However, due to 
a slight increase in background at large wave vectors (back- 
grounds tend to gradually, and reasonably isotropically, in- 
crease with higher Q due to multiphonon scattering), the 
slices [notably (a)-(c)] show excessive intensity at the high- 
est wave vectors. 
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Appendix A: Data corrected for Fe 2+ magnetic form 

factor 



To see more clearly the evolution of spin waves in 
Fei.osTe in Fig. 2 of the main text, we plot in Fig. SI 
constant energy slices of the data in Fig. 2 converted into 
units of magnetic response, \ (Q: w )- This is related to 
the measured cross-section by 



ki cPa 
VfdttdE 



TTg 2 ^L 2 B 1 — exp(— Tiuj/kT) 



(Al) 



where F(Q) is the Fe 2+ magnetic form factor, 
(7 r e) 2= 0-2905 b sr _1 , ki and kf are initial and final neu- 
tron wave vector and the g-factor was assumed to be 2. 
The exponential term is from the Bose factor, which can 
be neglected at the energies of interest, as all measure- 
ments were taken at 10 K. The left hand side of this 
equation is the raw data as plotted in Fig. 2. 



Appendix B: Heisenberg Hamiltonian 

For a given set of exchange energies (Ji , Ju, J<2, a , Jibi 
•hi Jz) as shown in Fig. 1(a), the dispersion of magnetic 
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excitations from the Heisenberg Hamiltonian with mag- 
netic order at (0.5,0,0.5) can be found by diagonalizing 
the Hamiltonian H [33] for every Miller index (h, k, I) 
value: 
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with 



" 2(Ji a + J 2a — Jib — Jib 

2J3[cos(2fc.j;) + cos(2fc a )] 



2J 3 

-2J Z 



+ Jz) 

cos k z ] 



C 



D 



+ e 
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k x = Tr(h + k) 
k y = Tr(h — k) 
k z = 2nl 



(B2) 



For a twinned sample, the full dispersion must be cal- 
culated at both (h, k, I) and (— k, h, I). The two Fe atoms 
per unit cell leads to two bands per wavevector, which af- 
ter twinning leads to a total of four bands. However, at 
least one of these bands has zero intensity at the wavevec- 
tor positions we compute the dispersion. 



Appendix C: The isotropy of J2 and the effect of Jj 
on the spin-wave fits 

As described in the paper, the dispersion data in 
the (h,0) and (1, As) directions were obtained from fit- 
ting wave vector peak positions of many constant en- 
ergy cuts. These two dispersion plots could then be si- 
multaneously fitted to the commensurate phase of the 
Heisenberg model to yield the exchange energies J\ ai Jib, 
Jia, J 2b, and J3. The parameters in the fits were con- 
strained to bounds such that magnetic order occurs at the 
(0.5, 0, 0.5) point, a requirement of the magnetic model 
(see the sections below and Ref. 26 for more details). 

To start with, the dispersion was fitted by varying all 
J ia, Ju, Jia, Jib, and J3 parameters. This yielded the 
fits in Fig. S2(a)-(b) which gave J u = -18.0 ± 7.4; 
Jib = -51.1 ± 3.4; J 2a = 22.6 ± 9.6; J 2b = 21.7 ± 3.5; 
J 3 = 6.77 ± 2.7 meV. This fit was robust, and as it 
was clear that J 2a and J 2o were very close, they were 
subsequently fixed to be equal and refitted. This fit 
(Jia = -17.5 ± 5.7; Ju = -51.0 ± 3.4; J 2 = J 2a = 
J 2b = 21.7 ± 3.5; J 3 = 6.8 ± 2.8 meV) [shown in Fig. 
4(a)-(b)] clearly yields parameters almost the same as 
the original fit parameters as expected, and is quoted as 
the main result in the paper. The low energy region of 



Fig. 4(a)-(b) is shown scaled up in Fig. S3 

Before discussing the necessity of the next-next-nearest 
neighbor term |34j . J3, we will first briefly compare our 
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FIG. S2: Dispersion extracted from data (solid symbols) and 
fits of the model dispersion (gray open circles with radius 
indicating intensity). Left panels show dispersion in (h,0) 
direction, and right panels show (l,fc) dispersion, (a)-(b) Fit 
performed with all parameters free, (c)-(d) Fit with fixed 
J3 = and J2 a — J2b- (e)-(f) Fit with fixed J3 = and 

J la = 0. 
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FIG. S3: Dispersion extracted from data (solid symbols) and 
model from best fit parameters (dotted line), zoomed in to 
emphasize low energy region of Fig. 4(a) and Fig. 4(b) respec- 
tively. Though for clarity the intensity of the model disper- 
sion is not shown in this plot, it can be seen from the marker 
sizes in Fig. 4(a)-(b) that the intensity of the branches to 
the right of 0.5 (and -0.5 respectively) in these figures would 
die away very quickly with increasing energy. Furthermore, 
at the lowest energies the two branches cannot be resolved. 
The outcome is that the branch to the right of the two figures 
are not observed in the data. 
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FIG. S4: Dispersion as shown by a slice of the raw data (mul- 
tiplied by energy transfer to enhance high energies) with a 
fixed h = 1, from the MAPS E t = 500 meV run. The disper- 
sion into (1,0) around 200 meV is clearly visible. 



dispersion with the theoretical exchange energies pre- 
dicted by Han et al. [50]. Figures S2(g)-(h) shows our 
experimental data compared with the dispersion calcu- 
lated for their theoretical values, and it is clear that the 
data and model are very dissimilar. The model can- 
not be reconciled with the data because in this theory 
J2b = —ha, whereas our data appears to be best de- 
scribed with Ji a = J2b- 

Following our successful fit with 32a — -j hb, it was 
important to check whether the small next-next-nearest 
neighbor parameter J 3 could be fixed to zero in order 
to further decrease the number of parameters in the fit, 
and see if this extra small exchange interaction is actu- 
ally necessary or not. The obvious starting place is to 
calculate the dispersion when using the above fit param- 
eters but with J 3 set to zero. However, the model disper- 
sion cannot be calculated in this case because this shifts 
the parameters out of the range allowed by (0.5,0,0.5) 
magnetism [26] . So instead, we performed another fit, 
whilst fixing J 3 to zero. In Figures S2(c)-(d), this fit 
can be seen, yielding J\ a — 6.7 ± 1.5; Ju = —31.7 ± 3; 
h = ha = hb = 9.6 ± 1.0 meV (and J 3 = 0). This 
fit describes the data below 100 meV well, but is not 
as successful at the high energies, underestimating the 
maximum energy of the band by around 50 meV. We in- 
terpret this fit as compensating the lack of an AF J 3 term 
by instead causing the J\ a term to be AF. The compro- 
mise is that Ji a cannot be too large or the high energy 
data cannot be described. 

One may finally ask whether the data could be fitted 
better if J 3 was fixed to zero and h a and hb could vary 
independently. In this case J\ a was found to stay small 
and so was fixed to zero for ease of fitting (results do not 
vary significantly with a fitted J\ a term). The fit is shown 
in Fig. S2(e)-(f) (with parameters J\b = —35.8 ± 2.9; 
J 2a = 14.0 ± 0.82; J 2b = 8.9 ± 1.3, and J Xa = 0; J 3 = 
meV). In this fit it can be seen that (i) h a and J2& are 
still quite similar, and (ii) the fit is also not as good as 
the original (J 3 > 0) fit, as it again underestimates the 
maximum energy. 
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FIG. S5: Instrument resolution convolved simulation of model 
with various exchange energy parameters, (a)-(h) fit with J3 
fixed to zero (and J2a = J2b), with dispersion shown in Fig. 
S2(c)-(d). (i)-(p) fit with J3 and J\ a fixed to zero, with 
dispersion shown in Fig. S2(e)-(f). (q)-(x) parameters from 
the Han et al. [20] prediction, the dispersion of which was 
shown in Fig. 4(c)-(d). For comparison with the simulation 
with best fit parameters shown in Fig. 5. J z does not vary 
the dispersion significantly (see main text of paper) and was 
therefore fixed for all plots to J z = 1 meV, the value chosen 
in Fig. 5. The intensity of the model in (q)-(x) was magnified 
by a factor of 1.5 in order to show the pattern more clearly. 



FIG. S6: Slices of the schematic phase diagram of the Ji — J2 — J3 model. The phases are defined in ref. 1261 among which 
the AF3 is the commensurate SDW phase observed in the experiment. ICA and ICB are incommensurate magnetic phases, 
and FM is a ferromagnetic phase. The red star in (a) shows the position in phase space of the best fit parameters shown in 
Fig. S2(a)-(b). The stars in (b) and (c) show the position in phase space of the other fits in Fig. S2(c)-(d) and Fig. S2(e)-(f) 
respectively. Note that the lower left corner gives the ferromagnetic phase, obtained from the (0.5, 0.5) order (AF1) phase by 
a sublattice symmetry S defined above. 



We conclude that the data is best fit with an isotropic 
J 2 (= J 2a = J2b)- The addition of a J3 exchange inter- 
action also seems to be important to describe the high- 
est energies whilst keeping the parameters in the regime 
with (0.5, 0, 0.5) magnetic order. For further insight into 
the similarity and dissimilarities between these different 
fits and models, we show instrument resolution convolved 
simulations with J3 fixed to zero [see Fig. S5(a)-(h) for 
simulation with parameters from the fit in Fig. S2(c)- 
(d), and see Fig. S5(i)-(p) for simulation with param- 
eters from the fit in Fig. S2(e)-(f)]. These alternative 
parameters show similar features to our data [compare 
with simulation of original ( J3 > 0) fit in Fig. 5], though 
do not extend high enough in energy. In Fig. S5(q)-(x) 
we have also included a simulation with the Han et al. 
(20] . theoretical parameters [see Fig. S2(g)-(h)] which 
have J 2 & = — J 2a- K is clear that these parameters do 
not describe the data well. 



Appendix D: Viewing the data from constant-Q 
perspective 

Performing energy cuts instead of (h, 0) and (1, k) cuts 
does not show the dispersion as clearly as Q-cuts because 
of the decrease in signal as energy increases. However 
the overall dispersion is clear when plotted as intensity 



multiplied by energy in constant-Q slices, see Fig. S4 It 
is still unclear in this slice whether the 250 meV almost 
non-dispersive band exists or not. 



Appendix E: The extended phase diagram 

In this section, we extend the phase diagram of the 
Ji — J 2 — J3 model (see ref. [2"rJl and Hamiltonian described 
in a section above) applied to FeTei-^Sen, system to the 
parameter region where Ji a ,ib can be negative (providing 
•ha > J2& > 0). First, it is only trivial to observe that 
the original model has the following two properties: 

U(S)^H(J la , Ji fc , J 2a , J 2b , J 3 )U(S) = 



H( — Jla,—Jlb,J2a,J2b,J3)> (El) 

U(R)*H{J la , Ji fc , J 2Q , J 2b , J 3 )U{R) = 

U(Jlb, Jla, J 2a , J 2 fc, J3), (E2) 



where S is the sublattice symmetry operation defined by 
S : S(i,j) — > (— Y +: >S(i, j) and R is a lattice rotation 
along the a-axis of angle tt, that is, R : S(i,j) — > S(j, i). 
This fact tells us we only need to focus in the param- 
eter region where Ji a > \Ju\ > 0, and the phases in 
other regions can be obtained by applying the appropri- 
ate symmetry operations on the ground state within this 
parameter region. For example, if one has Jij < Ji a < 0, 
first we obtain the ground state for the Hamiltonian 
H(—Jib,—Ji a ,...) and then transform the state under 
symmetry operation S x R. 

Second, as we have thoroughly discussed the phase di- 
agram with Ji a 16 > in the previous paper [26] . we only 
need to discuss the case where Ji n > while Jn, < in 
this note. This discussion is very simple because in this 
parameter region only one parameter, J 2o , is frustrated 
in the AFM3 phase, therefore the only possible instabil- 
ity is the one toward the ICB phase, the transition line 
of which is given by 



Jia — J 



16 



4J 



2b 



8J3. 



(E3) 



The full phase diagram is given in Fig. S6. 

We note that the effective exchange couplings are very 
close to the phase boundary between the (0.5, 0) AF or- 
dered state and two different incommensurate magnetic 
states with wave vectors (0.5 — S, 0) and (0.5 — 5, 0.5 + S). 
This is also consistent with the fact that the transi- 
tion from the (0.5, 0) AF state to the incommensurate 
(0.5 — 5, 0) magnetic state is observed in Fei +y Te by in- 
creasing y (refs. [HI IT5|) . There is also evidence that 
FeTei-^Se^ can have (0.5 — 5, 0.5 + S) incommensurate 
spin excitations [25l [35l Efi] . 
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